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Abstract 

We consider the "standard" spectrum of the active neutrinos (characterized 
by strong mass hierarchy and small mixing) with additional sterile neutrino, 
Ug. The sterile neutrino mixes strongly with the muon neutrino, so that 
i^n ^ oscillations solve the atmospheric neutrino problem. We show that 
the parametric enhancement of the <-> oscillations occurs for the high 
energy atmospheric neutrinos which cross the core of the Earth. This can 
be relevant for the anomaly observed by the MACRO experiment. Solar 
neutrinos are converted both to and Vg. The heaviest neutrino (~ v^) 
may compose the hot dark matter of the Universe. Phenomenology of this 
scenario is elaborated and crucial experimental signatures are identified. We 
also discuss properties of the underlying neutrino mass matrix. 







1. Introduction 

Reconstruction of tlie neutrino mass spectrum and lepton mixing is one of the funda- 
mental problems of particle physics. This problem may be solved with new generation of 
the neutrino experiments. 

It is believed that the most natural scenario of the neutrino mass and mixing is the one 
with strong mass hierarchy 

mg = (1 - 5) eV, ma = (2 - 3) ■ 10"^ eV, mi < ms (1) 

and small lepton mixing (comparable with mixing in the quark sector), so that z/g ~ z/i, 
~ 1^2, '^T ~ '^3- In this scenario the heaviest neutrino, u^, composes the hot dark matter 
(HDM) of the Universe, while the resonance conversion solves the solar neutrino 

(z/q) - problem [Q. This scenario (which can be called the "standard scenario") follows 
naturally from the see-saw mechanism with the neutrino Dirac mass matrix similar to 
the quark mass matrix: ~ m^p and the Majorana masses of the right handed neutrinos 
at the intermediate mass scale: 10^° — 10^^ GeV |Q. Short length z/g ^ and ^ ^'/z 
oscillations proceed via mixing of z/g and in the heaviest state z/3 The transition 
probability for the LSND experiment is small, P ^ 10"'^; It can be reconciled with data 
at 3(7 level only 0. 

In the standard scenario, there is no oscillation solution of the atmospheric neutrino 
{^atm) - problem [§]. Recently, the Super-Kamiokande and Soudan [|10| experiments have 



confirmed the existence of the problem. Moreover, observations of the zenith angle as well 
a,s L/E (distance/energy) dependencies of the muon neutrino deficit strongly indicate an 
oscillation solution of the problem. 

One can accommodate the oscillation solution of the Vatm- problem modifying the stan- 
dard scenario. A straightforward possibility is to diminish mg down to ~ (0.03 - 0.3) eV 
implied by the Vatm- data. For ~ (0.5 — 0.6) eV one can try also to explain the LSND 



result |TT| ignoring the zenith angle dependence in the atmospheric neutrino deficit. Further 
modification is needed to get sufficient amount of the hot dark matter in the Universe. In 
this connection the degenerate mass spectrum with mi ^ 1712 ~ ~ 1 — 2 eV has 
been suggested |T2|. Simultaneous explanation of the LSND result and the zenith angle 
dependence in the atmospheric neutrinos requires depart from the three neutrino scheme. 

In this paper we will consider a possibility to rescue the "standard scenario". The idea 
is to keep the standard structure with strong mass hierarchy (|l|) and small mixing for active 
neutrinos and explain the atmospheric neutrino problem by oscillations of z/^ to new neutrino 
state Us', ^ Us [l^. This new neutrino should be a sterile neutrino (singlet of SU2 x Ui) 
in order to satisfy the bound on the number of neutrino species from the Z°-decay width. 
Large — mixing can be related to a singlet character of Vg. 

The ^ Us oscillations as a solution of the atmospheric neutrino problem (see e.g. |ll3[ ) 



have been considered previously in the context of the model with mirror symmetry |14|. In 



contrast with [|14[ we introduce only one singlet state, and phenomenology of our scenario 



substantially differs from that in [|T1| . 

As we will show, in the suggested scenario the zenith angle dependence in the atmo- 
spheric neutrinos can be reconciled with LSND result at 3cr level (P ~ 10^'^). For relatively 
large mass of the new neutrino the predicted LSND probability can be enhanced up to the 
experimental value P ~ 3-10"^, however in this case no zenith angle dependence is expected. 

The paper is organized in the following way. In sect. 2 we describe the scenario. Then 
we consider predictions for the atmospheric (sect. 3) and solar (sect. 4) neutrinos, for 
the LSND experiment (sect. 5) and neutrinos from supernovas (sect. 6). Implications to 
the primordial nucleosynthesis are discussed in sect. 7. We comment on properties of the 
underlying mass matrix in sect. 8. 

2. Scenario 

The level and mixing scheme is shown in (fig. 1). The mass states z/i, z/2, z/3 have 
the hierarchy of the standard scenario (0). The mass of the additional state z/4 equals 



m4 ~ (0.5 — 7) • 10~^ eV, so that Am|4 ~ ml is in the range of solution of the Uatm- 
problem. The singlet neutrino, Ug, and the muon neutrino, z/^, mix strongly in the mass 
states 1^2 and z/4. We define the elements of the mixing matrix Uai which relate the flavor 
= (t'e, ^fi, Vs) and the mass z/j = (z/i, z/2, z/3, z/4) states as: 

T^a = , a = e,fi,T,s; i = 1,2,3,4. 

All mixings except for - Vg are small: ~ U^j_i ~ Us2 ~ Ugi = 0.5 — 0.8, and Ua,i <^ 1 
for others. 

3. Atmospheric Neutrinos 

A solution of the atmospheric neutrino problem is based, mainly, on z/^ ^ Ug oscillations. 
The muon neutrino state can be written as: 



- yl- U% - Ul^ ■ v' + ^7^3 • + ■ ^4 , (2) 

where v' is the combination of the light states v\ and z/2. According to our scenario the 
admixtures of the z/g in the heavy states z/3 and z/4 are small: f/e3, ?7e4 ^ 1- In fact, these 
matrix elements are strongly restricted by recent CHOOZ result [0. For small f/es and 
f/e4 the oscillation effects related to splitting between the states and z/2 are negligibly 
small. Note that splitting between these states in matter can be large enough, so that the 
oscillation length is smaller than the diameter of the earth. However, the effective mixing 
angle is strongly suppressed by matter. (The mixing may not be suppressed if [/e3 and ?7e4 
are large. In fact, f/e4 can be large, if < 3 ■ 10^^ eV. This possibility will be described 
elsewhere [0]). Thus, the task is reduced to oscillations in the (z/', z/3, z/4) system (see (0)) 
and the mixing elements [/^3, f/^4 are relevant only. Being in the eV - range, the state 
z/3 produces just the averaged oscillation effect, so that the z/^ survival probability can be 
written as 

P{y^ ^ y^) = (1 - It^^sH'A + If/^sl' • 

Here P2 is the survival probability for two neutrino oscillations in matter with parameters 
Am42 and 
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with the matter potential 

V = ±V2Gf 



liV^ + (i_fc)(iV,-iiV„ 



where k ~ |?^s4p/(|?^s4p + |^e4p) < 1- (The minus sign of V is for anti- neutrino.) Nn(t) 
and A^e(^) are the neutron and the electron number densities correspondingly. In the pure 
J^fi ^ t^s case one has k = 1, and V = ^j2GpNnl2^ so that the neutron density contributes 
only. With decrease of k the potential increases. 

If A'ml2 < 10^^ eV^, the oscillations are not averaged and P2 leads to the zenith angle 
dependence of the i/^ - deficit. For large the total effect is the combination of the 
averaged Ur oscillations and non-averaged u^j, Ug oscillations. 

The double ratio equals to survival probabihty: 

(/i/e)osc 
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as in the case of z/^ ^ u^. oscillation. In (H) (/x/ e)osc and (/i/ e) mc are the ratios of numbers of 
the /i-like to e-like events in presence of oscillations and without oscillations correspondingly. 
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z/o channels can follow from matter effect 



which influences the latter channel only. For maximal mixing the oscillation length and 
the effective mixing 6m which determines the depth of oscillations can be written as 
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For Am24 > 4 ■ 10"^ eV^ both for the sub-GeV and for multi-GeV (£' ~ 3 - 5 GeV) events 
we get ^ > 10 and the matter effect is negligibly small (< 10% in probability). Thus, in 
significant region of parameters (masses and mixing) at low energies the ^ Vg oscillations 



reproduce results of the z/^ ^ z/^ oscillations, as far as the charged current (CC) interactions 
are concerned. (The number of the tau-lepton events is small.) 

The matter effect can be important, for multi-GeV events, if Am24 < 3 ■ 10"'^ eV^. In 
particular, the zenith angle dependence of ratios can be modified. Since matter suppresses 
the z/^ ^ Ug oscillations, one expects weaker oscillation effect for the upward-going events, 
and therefore, flattening of the zenith angle, G, dependence. In particular, R{cosQ = 
— 1) > 1/2 in the z/^ ^ z/g case, whereas for the z/^ ^ z/r case one gets i?(cos = —1) = 1/2 
provided the mixing is maximal. 

The matter effect becomes important for through-going and stopping muons [|l^ pro- 
duced by high energy neutrinos even for large Am24. With increase of energy (decrease 
of Am^) ^ — > and the oscillation length increases approaching the asymptotic value de- 
termined by the potential only: ^ 27i/V. In the same time, the effective mixing angle 
decreases, so that the oscillation effects become weaker. 

The zenith angle dependence of the survival probability for fixed energies is shown in 
fig. 2. According to fig. 2 the dependence has two dips: the wide dip with minimum at 

cos© ~ (—0.5 0.4) and the narrow one with minimum at cos0 ~ —0.9. With increase 

of energy the depths of dips (peaks in the transition probability) decreases, whereas the 
shape of dips, and in particular, positions of minima change weakly. This behavior can be 
understood in the following way. The positions of minima and maxima are determined by 
the phase of oscillations: 



For E > 20 GeV and Am^ < 5 ■ 10"^ eV^ we have ^ < 0.2, so that the length of oscillations 
and therefore the phase are determined by the potential V and only weakly depend on energy. 
As a consequence, the phase of oscillations is fixed by the zenith angle. According to the 
model of the Earth ||T^, the trajectories touch the core at cos 9 ~ —0.8. For cos 6 > —0.8 
neutrinos cross the mantle only, and the wide dip at cos0 > —0.8 is due to the oscillations 
in the mantle. At cos O ~ —0.4 the phase $ equals vr which corresponds to the minimum of 
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p. It turns out that at cos0 ~ —0.8 the phase is $ = 27r and the oscillation effect is zero. 

At cos O < —0.8 neutrinos cross both the mantle and the core. That is, the narrow dip is 
due to oscillations of neutrinos whose trajectories cross the core. In this dip the oscillation 
effect is even stronger than in the wide dip, in spite of larger density of the core. This dip (or 
peak in the transition probability) is the parametric resonance peak [|TB],|T^ : the parametric 



enhancement of the z/^ ^ Us oscillations occurs for neutrinos which cross the core of the 
Earth. Indeed, these neutrinos cross three layers with approximately constant densities: the 
mantle, the core, and again the mantle, divided by sharp change of density on the borders. 
Let us denote the phases acquired in these layers by $mi, ^c, ^m2 (obviously, = ^m2)- 
It turns out that for cos© ~ —0.88 the following equality is realized: 

$ml ^ $c ~ ~ vr . (7) 

That is, for each layer the size of the layer coincides with half of the averaged oscillation 



length: /mj/2 ~ L^. This is nothing but the condition of the parametric resonance ||T8|,|T9 
The amplitude of oscillations is enhanced. Maximal transition probability which can be 
achieved in the case of three layers equals: 

Pi, = sin2(4C""*-2C'''), (8) 

where ^J^**"* and 6'^'^'^ are the averaged mixing angles in the mantle and in the core corre- 
spondingly. Let us stress that for sufficiently large energies the equality (0) does not depend 
on neutrino masses. The equality is determined basically by the density distribution in the 
Earth and by the potential which in turn is fixed by the channel of oscillations and the 
Standard Model interactions. The equality is fulfilled for oscillations into sterile neutrinos 
only. For the z/^ <-> z/g channel the potential is two times larger and (|^ is failed. Note that 
for cosO = — 1 the oscillation effect is small. 

Qualitatively the zenith angle dependence shown in fig. 2 is similar to that observed 
in the MACRO experiment ||2^. There are two dips with minima at cos0 = (—0.4 — 



-0.6) and cos© = (— LO 0.8). In our interpretation the latter dip could be due to 
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the parametric enhancement of the z/^ ^ Vg oscillations of neutrinos which cross both the 
mantle and the core of the Earth. Similar features, although not so profound, have been 



observed in the Baksan experiment ||2T| which has energies of the detected neutrinos similar 
to those in the MACRO experiment. In the Super-Kamiokande experiment the effect is 
expected to be weaker, since the energy threshold of the through-going muon detection and 
corresponding energies of neutrinos are higher Detailed comparison of the effect with 



data will be given elsewhere 



Another consequence of the z/^ ^ Vs oscillations is that the region excluded by the 1MB 
analysis of the stopping and through-going muons is shifted by factor ~ 3 to larger values 
of Am? |13|. This allows one to reconcile the excluded 1MB region with the preferable 
Super-Kamiokande domain. 

Crucial check of the ^ Vg oscillations can be obtained from studies of events produced 
by the neutral current {NC-) interactions of the atmospheric neutrinos. The number of 
the iVC-events is suppressed by the z/^ Ug oscillations and unchanged by the ^ 
Vr oscillations. It is difficult to detect the elastic neutral current interactions. In this 
connection, it was proposed to study tt^'s produced, mainly, by the neutral currents pB| , pi 



m 



P3| the long base line experiment K2K was discussed). Subsequent decay n — > 77 
is identified as the two showering event with certain invariant mass. Practically, one can 
measure the number of vr^-events, N^o, and find the ratios N^o/N^ and N^^o/N^, where A^e 
and A"^ are the numbers of e-like and /i-like events induced, mainly, by the charged currents. 
This eliminates uncertainties related to absolute values of the atmospheric neutrino fluxes. 
Another proposal |^ is to measure the ratio of A^^o and the number of events with 1 vr- 
production by the charged currents. The latter dominate in the two-prong (ring) events 
sample, N2ring, and experimentally it is possible to determine the ratio 



This allows one to diminish uncertainties related to cross-sections. Also one can study 
the ratio N^o /N mum-ring- Since ^ oscillations suppress equally the CC- and NC- 
interactions of the muon neutrinos, while the z/^ Vr oscillations suppress the CC- events 
only, one expects 

Similar inequality is expected for N^o/Ng. For the survival probability P ~ 0.6 the difference 
can be as big as (30 — 50)% Thus with about 300 7r°-events which will be detected 



by the Super-Kamiokande experiment within 2-3 years one will be able to disentangle 
solutions. Additional check will be done by the long baseline experiment K2K (KEK to 
Super Kamiokande) In particular, one will be able to compare the number of tt^'s 

produced in the front (close to neutrino production target) detector and in the Super- 
Kamiokande detector. However, the number of expected in the SuperKamiokande is 
expected to be small ( ~ 30 ). 

4. Solar Neutrinos 

Due to strong mass hierarchy, a dynamics of conversion of the solar neutrinos is reduced 
to one Am^ task. Indeed, the electron neutrino state can be written as 

i^e = cos (j) ■ u'^ + sincj) ■ , (10) 

where 



sin0^V^f/23 + t/2, , (11) 
and is the combination of the light states: 

= cos a; ■ z^i + sina; ■ 1^2 (12) 

with 

cosu = I ^'^^ = . (13) 
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The uh is a combination of the heavy states: 

= 7^ + • . (14) 

For evolution of the solar neutrinos produced as v^, we get the following picture. The states 
z/3 and z/4 decouple from the system, thus leading to the averaged oscillations result. The 
state z/g converts resonantly to its orthogonal state: 

= coscj • z/2 — sincj ■ Ui . (15) 

Using this picture we find that the survival probability for the solar Uf. can be written as 

Piu, i^e) = cos" (I) P2MSW + \Ue4\^ + Pesl^ ■ (16) 

(Here we also take into account that oscillation effects in z/3 — 1/4 system are averaged out.) 
P2MSW is the z/g survival probability of z/^ — > z/^. resonance conversion. The amplitude of the 
probability can be obtained by solving the evolution equation 

for the two neutrino system u = (z/g, I'x)'^- The effective Hamiltonian 
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can be found from the whole 4- neutrino Hamiltonian. In ( |TB|) Am2i = m2 — mf, and is 
determined by 
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^^^'^= I.. 2" 12 • (19) 



Note that for small ?7^3 the angle ip coincides practically with Oatm (i|) responsible for oscil- 
lations of the atmospheric neutrinos. In ([18|) the effective density equals 

N'ff = Ne cos^ - ]^Nn{sm^ ^-\Ue4\^ cos^ V^)- (20) 



If mixing of Ve in heavy states is negligibly small, the whole survival probability is reduced 
to P2MSW according to (|I6]). 



In our scenario new features of the matter effect appear in comparison with pure active (z/g 
Vfj) or pure sterile (z/g — >■ z/^) neutrino conversions. In particular, for non-zero admixture 
of the z/g - flavor in the state z/4 {Uei 7^ 0) the off-diagonal elements of the Hamiltonian (p!8D 
depend on matter density, so that the mixing can be induced, at least partly, by the matter 
effect. Let us consider the ratio of the vacuum and the matter contributions to mixing: 

f/e4 sin 2i) ■ y^GpK^' 



Rm/v 



Am 



(21) 



21 



E 



sin 20; 



Here the neutron density is taken in the resonance = N^'^^ which is justified for small lo 
(narrow resonance). Determining A^^^* from the resonance condition {H22 = 0, where H22 is 



the 22-element of the Hamiltonian ([T8|) ), we find 



R, 



■m/v 



e4 



4 tan 2u 



COS^ — 2 (^^^^ ^ ~ \Uei\^ COS^ 1p) 



u, 



e4 



4 tan 2uj 



(22) 



where the last equality holds for small and ip ~ tt/A. Here r^j = Nf,/Nn is the ratio of 
the electron and neutron densities in resonance. The ratio increases from 2 in the center of 



the Sun to about 6 at the surface |25]. Using the upper bound f/e4 < 0.2 from the CHOOZ 



experiment |15] and the lower bound, tan2u; > 0.06, implied by a small mixing solution of 
the z/Q-problem we get from (^) Rm/v ^ 0.4. That is, the matter effect is smaller than 
the vacuum effect, although the difference is not large. The matter effect can both enhance 
and suppress the vacuum effect, thus enlarging a range of possible vacuum mixing angles: 
sin^ 2u = (0.5 — 3) sin^ 2uq; here c^o is the angle without off-diagonal matter effect. 

If the original boron neutrino flux is smaller than the SSM one, then smaller total mixing 
is needed and Rm/v can be of the order one. In this case the matter effect can give dominant 
contribution or even substitute the vacuum mixing in ([l^). 
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For Ue4 <C 0.2 the matter term in the off-diagonal element of H2 is much smaller than the 
vacuum term, and P2MSW reduces to usual 2z/ survival probability characterized by u, Amg^ 
and N'^-l'^ ^ cos^ (pN^ — \ sin^ ipNn (see (pOD). Using the reactor bound on the admixture of 
Ve in the state z/3: cos^ > 0.95, we get Ar^*^"^^ ^ ^eff ^ ^^^^^ ^sterile = - Nn/2 
is the effective density for pure Ve ~^ conversion. As a consequence, in our case the 
survival probability at the adiabatic edge of the MSW suppression pit (survival probability 
as a function of E) is in between the probabilities for the 2v sterile and 2v active cases; the 
non-adiabatic edges (determined by the derivative d[ln{N'^^^)]/ dx) practically coincide for 
all three cases. 



Main feature of solution of the solar neutrino problem in our scenario which can be used 
for its identification is that solar neutrino flux should contain both and neutrinos. 
Moreover, due to large - Vg mixing in the z/2-state implied by the atmospheric neutrino 
data the transition probabilities P(z/e — > Vg) and P{ve ^n) should be comparable. (The 
z/g — > v-r probability is small, since v-r is weakly mixed with other neutrinos.) We find: 

-P(z^e ^ ^'m) ~ COS^ Sin^ Ip Sin^ (t)P2MSW + cos^ ^"(1 - P2Msw) 

+ sin sin 2ip tan 2c<j(i — P2msw) + sin^ sin^ ifj , 



Pfz/p, — > z/.J ^ cos 



k2 



cos^ ip sin^ (j)P2MSW , , 

L (24) 

+ sin^ -0(1 — P2Msw) — sin sin 2-0 tan 2ijj{\ — P2MSW) + sin^ cos^ ^ . 
In the limit of small uj and small we get from (^, |2^) : 

P{ye Vy) ~ COS^ V^(l - P2MSW) . . 

(25) 

That is, the relative contributions of these channels are determined by the angle ■0. In 
particular, for sin^ -0 ^ 0.5 we have P{ue u^) ^ P{i'e ^s)- Taking sin^ 2-0 ^ 0.70 — 1.0 
as it follows from the atmospheric neutrino data we find 
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This ratio does not depend on the neutrino energy (fig. 3). 

In what follows we will refer to this solution as to the mixed z/g — > z/^, z/^ solution. As we 
will see, its properties are intermediate between properties of the pure fiavor z/g — > z/^ and 
pure sterile Ve solutions. 

Some general tests of presence of the sterile neutrinos in the solar neutrino fiux have 
been elaborated in Transition to sterile neutrinos suppresses the number of neutral 

current events. Let us consider a modification of the charged-to-neutral current events 
ratio, (CC/NC), which will be measured in the SNO experiment. Let us introduce the 
double ratio: 

/CC\ /CC\ 
[ncJosc ' KNCJo ' 

where {CC/NC)o and [CC /NC)osc are the ratios in absence of oscillations and with oscil- 
lations. For three possible solutions of the z/Q-problem we get: 



<Pee> 
< Pee > 
1- < Pes > 



ye 



Ve 



(27) 



where < ... > means averaging over energy spectrum of the boron neutrinos. According to 



rd{l^e l^fj,) < rd{ve v^,, Us) < rd{Ve Vs) ■ 



(28) 



For small mixing solution and the SSM boron neutrino fiux [25|, we find intervals 





0.3 - 0.4 


^e - 








0.5 - 0.8 


Z^e - 


^ Z/^, Z^s 


(29) 




1.07 - 1.09 


l^e - 


z/s 





which are well separated from each other. (In (^) difference of energy dependences of the 
CC- and NC- current cross-sections is taken into account, which leads to > 1 in the last 
line P3].) Thus, if the original flux of the boron neutrinos, Fb, is known, it will be easy to 
disentangle solutions by measuring r^. However, uncertainties in Fb make the task to be 
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ambiguous. If, e.g., Fb is smaller than the one in the SSM Fb ~ 0.7F§^^, then for 
the (i^e ^fi) solution we get ~ 0.4 — 0.6 which overlaps with the interval expected for 



Due to presence of both and in the final state, the effect of the A^C-interactions in 
the z/e— scattering (the Super-Kamiokande experiment) is intermediate between the effects 
in the pure active and pure sterile neutrino cases. Correspondingly, the sin^ 2u! — Am2i 



region of solutions is intermediate between regions for pure active and sterile ||26| neutrino 
conversions. This gives sin^ 2uj = (4.5 - 11) ■ 10"^ and Am^^ = (3 - 9) ■ 10"^ eV^ in the SSM 
framework. 



Let us consider the distortion of the energy spectrum of the recoil electrons. In 
the SNO detector (i/grf epp) the distortion is determined by the probability P(z/e — >• 
z/e)~ cos'^ (j)- P2MSW- The factor cos'* (p leads to smoothing of energy dependence of the prob- 
ability, and therefore the distortion of the recoil electron spectrum can be weaker than in the 
pure Ue ^fj, or z/g — > case. However for sin^ 2(f) < 0.2 this effect is small (cos* cj) ~ 0.95). 
The Super-Kamiokande experiment (z/e ve) is sensitive to both CC- and NC- inter- 
actions. The A^C-scattering of leads to smoothing of the spectrum distortion. The 
transition to sterile neutrinos removes this effect. Therefore in the case of mixed conversion, 

— ^ v^iii-'si the spectrum is distorted stronger than in the Ve i^fj. case but weaker than 
due to z/e —>■ Vs conversion [^] (fig. 4), although the difference is small. 



Let us discuss the day-night effect |^0|. The z/g — > z/s,z/^ gives weaker asymmetry than 
i^e i^^i, but stronger asymmetry than z/g — > z/^. The last two cases were studied in and 



An important signature of our scenario is the day-night effect for events induced by the 
neutral currents, A^^*". (This in pure z/g z/g conversion case has been studied in |2^.) Let 
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us introduce the day-night asymmetry as 



A 



NC _ ^^Af ^"D 



(Similarly, the asymmetry for CC-events, A^'Jj^, can be defined). The asymmetry can be 
used to distinguish our scenario from pure 2i/-cases. Indeed, A^'j^j^ is zero for the 



z/p z/„ 



conversion and non-zero in presence of sterile neutrinos. The intermediate case z/^, z/^) 

can be distinguished from the pure sterile case by comparison of asymmetries in the charged 
currents, A'j^j^, and in the neutral currents. For pure sterile case: A'^^j^ — A^'j^j^ p6 



In the (z/g — > z/^, Ug) case we predict A^Jj^ — A^j^j^ > 0. Studies of these asymmetries will 
be possible in the SNO experiment. 

5. Short range oscillation experiments. LSND. 

In our scenario, two neutrino states with large masses, z/3 and z/4, can be relevant for 
short range experiments. The probability of z/^ ^ Ug (z/^ ^ z/g) oscillations can be written 



as 



P{U, Z/,) ^ 4|f/,3nf/e3psin2 \^]+ W\Ue4\' siu' ^ + 



ARe[U:,U^3U,^U*J ■ sin 



(m| — ml)t ^ 6 



4E 



- sm^ + - - sm^ — 

V 2 \ iE 2 



(30) 

where 5 = arg{U*^Uf^3Ue4U*^). For real mixing matrix we get 5 = 0. The first two terms 



in ([30| ) correspond to "indirect" oscillations related to levels z/3 and z/4 correspondingly. 
The third term is the result of interference of oscillations due to and u^. For ~ 0(eV) 
we get the following results depending on value of the mass 7714: 

1) The mass < 0.3 eV is too small to contribute to the LSND effect and the problem 
is reduced to the one level problem with probability described by the first term in (0). For 
values of and Ues similar to those in quark sector we get negligibly small probability for 
LSND: P < 10^^. If f/^3 and f/gs are at the upper experimental bounds, then P < 10^^ 
which could be reconciled with data at 3a-\eve\ For < 0.1 eV, one can reproduce a 
desired zenith angle dependence in atmospheric neutrinos. 
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2) For 777,4 > 0.3 eV both and 1/4 can give sizable effect to tfie sliort range oscillations. 
Also interference between two different modes of oscillations becomes important. In certain 
ranges of parameters the interference is constructive, thus leading to substantial enhance- 
ment of the probability. For instance, iiml — 1 eV^ and m| = 0.5 eV^, the probability P{i^ij, 

Ug) is enhanced by a factor of more than 4, as compared to the one level case. Now the 
probability P ~ 2.2 • 10~^ turns out to be in a good agreement with LSND result. In this 
case one also predicts <-> i/^ and <-> oscillations at the level of present experimental 
bounds. However, no zenith angle dependence of atmospheric neutrino deficit is expected. 

6. Supernova neutrinos 

Since i/^ and Us arc (almost) maximally mixed, their level crossing occurs at zero (small) 
density. Let us introduce the eigenstates of system (i/^, Ug) in medium u!^^, 1/'^^, in absence 
of mixing with i/e and i/r- At zero density u!^^ ^ 1/2 and 1/4^ ~ 1/4. Using the eigenvalues of 
1/2^ and it is easy to find that whole system {ug and Ur included) has four resonances 
(level crossings) in the neutrino channels (fig. 5): (i) {ut - u'^^ crossing at density (h) 
(i/g - fr) crossing at density A^er, (hi) {^e - J^im) crossing at A^e4, and (iv) - 1^2^) crossing 
at Ne2- The mass hierarchy leads to hierarchy of the resonance densities (see fig. 5): 

Nr4 > Ner > Ne4 > Ne2. 

There is no resonances in the antineutrino channels. For ma < 5 eV we get niNNrA < 
10^ g/cm^ (here mjv is the nucleon mass). That is, the Vr i^s resonance conversion 
occurs outside the core. Moreover, the effective v^. — Ug mixing in the core is very small: 
sin^ 2^^ Ri 4:\Us'i\^ -NrA/Ncore , and generation of the sterile neutrinos is strongly suppressed. 
Consequently, there is no influence of the sterile neutrinos on the gravitational collapse and 
coohng of the core. 

The nucleosynthesis of heavy elements due to the r-pro cesses in the inner parts of su- 
pernovas implies that Vr — > conversion does not occur efficiently above certain densities 
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3^ , so that there is no electron neutrinos with hard spectrum in the range (10^ — 10^°) 
g/cm^. The mass spectrum under consideration reahses the possibihty suggested in p 



which allows one to enhance the r-processes. Indeed, suppose that transitions v-r Vg and 
z/g Vr due to both high density resonances are efficient. In this case in the region of the 
r-processes the v^, - flux is absent or strongly suppressed (see fig. 5), thus creating neutron 
rich medium needed for production of the heavy elements. 



Let us consider consequences of transitions for neutrino fluxes detectable at the earth. 
Suppose ffist that all resonances are effective. In this case the following transitions occur: 

(31) 

Thus one expects (i) that z/g-flux will have hard spectrum corresponding to the original z/^ 
- spectrum; (ii) in contrast, z/^'s will have soft spectrum (corresponding to the original z/g- 
spectrum); (iii) the ffiix of z/^'s is produced, of the order of ~ 1/12 of the total flux. 

Suppose now that low density resonances u^. — z/2m and z/g — z^4m are unefficient (which 
may occur in compact stars with small mass and large density gradient in the outer layers). 
Then transitions proceed as follows: 

Ur ; l^^i-^ V2 (^ Vp,, l^s)] Z/^ ^ Z/g Z/4 U^, Ug). (32) 

Thus z/g's disappears which can be established in experiment, z/^'s will have soft spectrum 
of the original z/g, the z/g-fiux will compose ~ 1/6 of the total flux. There is no transitions 
in the antineutrino channels. 

Finally, suppose that Ug —>■ Vr resonance is not efficient because of smallness of t/gs. The 
following transitions take place: 

z/g z/4 (^ z/^, Vs) ; z/^ z/g; ^2 (~ v^, v,). (33) 

In this case the z/g-fiux in the r-processes region is unchanged. At the detectors, however, the 
the z/g-fiux will have hard spectrum. The z/^- and v^- fluxes will have both a hard component 
due to the z/^ convertion and a soft component due to the z/g z/4 convertion. 
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7. Primordial Nucleosynthesis 

The oscillations i^fi ^ i^s, i^fj. ^ with parameters implied by the atmospheric neutrino 
problem result in appearance of the equilibrium concentration of sterile neutrinos in the 
early Universe. Therefore, during the time of nucleosynthesis, t > Is, the effective number 
of neutrino species is N,^ = 4. This disagrees with bound < 2.6 based on low values 



of deuterium abundance |24[. Four neutrino species are, however, in agreement with recent 
conservative estimations N^, < 4.5 [25|,0 which use large abundance of '^He. In view of 
these uncertainties let us comment on possibility to suppress generation of sterile neutrinos. 
It was shown that leptonic asymmetry in active neutrinos Lj, = {n„ — np)/n^ > 7 ■ 10~^ at 
the time t > lO^^s is enough to suppress ^ oscillations. Such a produces 
the potential due to u — u scattering: V ~ GpL^n^ which suppresses the effective mixing 
angle in matter: sin^ 26'^ ~ sin^ 2^Am^/Tl^ (T is the temperature). Although there is 
no strong restrictions on leptonic asymmetry one would expect that L,, is of the order of 



the baryon asymmetry: L,^ ~ 10 It is argued |^ (see also |^) that large can be 



produced in oscillations of the tau neutrinos to sterile neutrinos: i^s, i^r in earlier 

epoch (t < 10^^ s). In presence of small original leptonic asymmetry, Ly ~ 10^^°, these 
oscillations can lead to exponentially fast increase of (the task is non-linear). For m-^ > 1 
eV, the increase occurs at T ~ 16 MeV. The lower value of mg required by this mechanism 
is approximately proportional to Aml2- Thus, for Aml2 = 3 ■ 10^^ eV^ we get from ||3^ 
m| > 6 eV^ and sin^ 29rs < 8 ■ 10~^. For Am^g = 10"'^ eV^ the corresponding bounds are 
mg > 2 eV^ and sin^ 26't-s < 4 • 10^^. For larger angles the u^- oscillations themselves 

produce the equilibrium concentration of sterile neutrinos. These conditions can be satisfied 
in our scenario. 
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8. Properties of mass matrix 

Let us find the elements of the neutrino mass matrix in the fiavor basis, ||mQ/3||, {a, (3 = 
e,fi,T,s) which lead to the scenario under discussion (fig. 1). Suppose first that mixing of 

is very small and it does not infiuence masses and mixing of lighter neutrinos. Then in 
the rest system: (z/^, z/g, z/^) the sub-system (z/s,z/^) with bigger masses and approximately 
maximal mixing dominates. This means that we can find masses and mixing of and z/^ 

(2) 

considering the 2x2 sub-matrix | |, with a, (3 = s, n only. The effect of mixing with z/g 
gives small corrections. 

Main feature of our scenario is the maximal mixing of two states with at least one order of 
magnitude mass hierarchy: e = 1112/111^ = 0.03 — 0.10. Maximal mixing and strong hierarchy 
imply that sub-matrix has small determinant and all its elements are of the same 

order. Indeed, we find the following relations: 

m^s tan2^(jtm 

m^^niss - _ e 

{m^^ + rUssY (1 + e)2 

and m4 m^^ + mss- For maximal mixing the masses are determined by 

™.^»..™..i±^. (34) 

For 0.03 — 0.1, the eq. (H) gives m^s ~ (1.05 — 1.2) ■ m^^ or "rriss ~ —(0.95 — 0.8) ■ m^g. In 
the case of non-maximal mixing larger spread of masses is possible. For sin^ ^Oatm = 0.9 we 
find: niss = fri^^ + ^^is, and m^^ : m^s '■ ^ss = 0.67 : 1 : 1.3. Thus, a spread of the 
elements can be within factor of 2. Taking m4 ^ 2m^s, we get: m^s ~ (1 — 3) ■ 10^^ eV. 

Diagonalizing the sub-matrix m^^^ one can find the mass matrix for the light neutrino 
system (z/e, ^'2)1 where z/2 is the neutrino with mass m2- From this we get the mixing of the 
light neutrinos: 

sin6'0 ^ —{rrie^cosOatm - mesSmOatni) ^ (3 - 5) ■ 10"^ , (35) 
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and if there is no strong cancellation in (pS]): 

T^e/i, T^es < sin ^0 -1712 ~ 2 ■ 10~^ eV . 

The structure of the mass matrix can be substantially different if the mixing of the tau 
neutrino is not small. In particular, the desired mass rriss can be generated by this mixing: 
trLss ~ —{msT-Y/m^T-, if, e.g., m^,- ~ 0.1 eV and 171^^. ^ (1 — 2) eV. This corresponds to 
rather large Ur — mixing parameter: sin^ 20 rs ~ (2 — 4) ■ 10^^ for which the equilibrium 
concentration of Vs is produced in the early Universe already in ^ Ug oscillations. The 
mechanism of generation of the lepton asymmetry |^ does not work. 

As follows from the above estimations of the mass terms, the mass matrix for the active 
neutrinos can have usual hierarchical structure with small mixing. It can be easily generated 
by the see-saw mechanism with linear hierarchy of masses of the right handed neutrinos. 
Also Ug can have the hierarchy of couplings with different generations: rries <^ ttLus -C rrirs. 
Some elements can be zero: e.g., nies = = 0. 

9. Conclusion 

We have considered minimal modification of the "standard" scenario for the neutrino 
mass introducing one additional sterile neutrino. The sterile neutrino mixes strongly with 
muon neutrino, so that the ^ Ug oscillations solve the atmospheric neutrino problem. 

We show that the parametric enhancement of the z/^ ^ oscillations occurs when 
high energy atmospheric neutrinos cross the core of the Earth. This effect can be relevant 
for explanation of the anomaly in the zenith angle distribution observed by the MACRO 
experiment. Suppression of the A^C-events in the atmospheric neutrinos (tt^ - like events) 
is another signature of the scenario. 

The scenario can supply the hot dark matter of the Universe. 

The solar neutrino problem is solved by z/g conversion to z/^ and Ug . The solution has char- 
acteristics (spectra distortion, day-night effect, NC/CC ratio) being intermediate between 
characteristics of — t'^ and z/g — z/^ conversions. 
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The probability of — > Pg oscillations in the LSND experiment can reach O(10~^). It 
can be further enhanced if new neutrino state has the mass 777.4 > 0.3 eV. (In this case, 
however, the atmospheric neutrino deficit will have no zenith angle dependence.) 

Disappearance of the i/g-flux from supernovas is another possible signature. 

The scenario can be well identified or rejected by the neutrino experiments of new gen- 
eration (Super-Kamiokande, SNO, Borexino, ICARUS and others ). 
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Figure Captions 



Fig. 1. Qualitative pattern of the neutrino masses and mixing. Boxes correspond to differ- 
ent mass eigenstates. The sizes of different regions in the boxes determine flavors of mass 
eigenstates, If/^jp- White regions correspond to the sterile flavor, hght and dense shadowed 
regions fix the electron and muon flavors correspondingly; admixtures of the tau flavor are 
shown in black. 

Fig. 2. The zenith angle dependence of the survival probability P{i'ij, — i^n) for atmospheric 
neutrinos for different values of the neutrino energy (figures at the curves in GeV) and 
Am^ — 5 • 10~^ eV^. Solid hues correspond to <-> Vg oscillations, dashed line is for 
^ oscillations. 

Fig. 3. The conversion probabilities of solar neutrinos as functions of E/Am? for 
sin^ 2u! = 0.01, cos^ ip = 0.5 and sin0 = —0.2; P{i'e ~^ ^e) is shown by solid line, Pi^e ~^ ^n) 
- dotted hne, P{ve ^^s) ~ dashed hne. 

Fig. 4. The expected distortion of the recoil electron energy spectrum in the Super- 
Kamiokande experiment. The ratio of numbers of events with and without conversion 
is shown by histograms: bold solid line corresponds to i^e — > i^s, i'^ conversion considered in 
this paper; solid line is for u^. — >■ conversion, and dotted line is for z/g i^^ conversion. 
The following values of parameters were used: Am^ = 5 • 10~^ eV^, sin^ 2uj — 8.8 • 10~^, 
and (for our scenario) cos^ V — 0-5 and sin = 0. Also the Super-Kamiokande experimental 
points from 306 days of observation are shown (statistical errors only). 

Fig. 5. The level crossing scheme for supernova neutrinos. Solid lines show the eigenvalues 

of i^3mi ^Am, ^2m and Pim- Black dots represent 4 resonances. Dashed lines correspond to 
energies of i/g, (straight lines) and i'2m- 
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